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Using coupled Ginzburg-Landau equations, the dynamics of hexagonal patterns with broken chiral 
symmetry are investigated, as they appear in rotating non-Boussinesq or surface-tension-driven con- 
vection. We find that close to the secondary Hopf bifurcation to oscillating hexagons the dynamics 
are well described by a single complex Ginzburg-Landau equation (CGLE) coupled to the phases of 
the hexagonal pattern. At the bandcenter these equations reduce to the usual CGLE and the system 
exhibits defect chaos. Away from the bandcenter a transition to a frozen vortex state is found. 



The interplay between theory and experiments is at 
the heart of the progress made in recent years in the 
study of pattern-forming systems far from equilibrium. 
Thus, investigations of specific physical systems that al- 
low a quantitative comparison with precise experiments 
are essential to check available theories and often also to 
suggest new theoretical venues. This has been particu- 
larly successful for ordered patterns (e.g. in Rayleigh- 
Benard convection and Taylor-Couette flow ||J). For 
disordered patterns and spatio-temporal chaos a num- 
ber of different systems have been under investigation 
with the goal to achieve detailed comparison between 
experiments and theoretical work. Spiral-defect chaos 
in Rayleigh-Benard convection with low Prandtl number 
exhibits very rich dynamics (e.g. Q). However, they do 
not arise at small amplitudes and are therefore not ac- 
cessible by weakly nonlinear theory. In the presence of 
rotation the Kiippers-Lortz instability of convection rolls 
induces another type of spatio-temporal chaos (e.g. ||). 
Although it occurs directly at onset a systematic theo- 
retical treatment within weakly nonlinear theory has not 
been successful due to the isotropy of the system. In 
anisotropic electroconvection of nematic liquid crystals 
spatio-temporal chaos jl| can be systematically described 
by coupled Ginzburg-Landau equations. The derivation 
of these equations from the microscopic equations is, 
however, extremely involved and not yet complete [^). 

From a theoretical point of view the most attractive 
and therefore most extensively studied canonical equa- 
tion exhibiting spatio-temporal chaos is the complex 
Ginzburg-Landau equation (CGLE). It describes the on- 
set of oscillations in a spatially extended system. De- 
spite its simplicity it exhibits an extraordinary variety of 
complex dynamics, including phase chaos, defect chaos, 
and frozen structures JoT-pl . In the one-dimensional case 
many aspects of its dynamics have been observed exper- 
imentally (e.g. localized solutions |15|). In two dimen- 
sions, however, no detailed comparison of any of the the- 
oretical regimes of complex dynamics with experiments 
is available. 

In this Letter we suggest that weakly nonlinear 
hexagon patterns in rotating convection are a good can- 
didate to compare the theoretical results for spatio- 



temporal chaos in the two-dimensional CGLE with ex- 
periments. Due to the rotation the hexagons typically 
undergo a transition to oscillating hexagons, which are 
described by a complex Ginzburg-Landau equation cou- 
pled to two phase equations. We show that in a suf- 
ficiently large system the oscillating hexagons exhibit a 
state of defect chaos that is well described by the usual 
two-dimensional CGLE. Depending on the wavenumber 
of the underlying hexagons we find a transition to a 
frozen vortex state, as it is also the case in the CGLE. 

We consider three coupled amplitude equations, de- 
scribing the dynamics of hexagonal patterns with bro- 
ken chiral symmetry, as they appear in rotating non- 
Boussinesq or surface-tension-driven convection. These 
equations can be obtained from the corresponding phys- 
ical equations (e.g. Navier-Stokes) for small amplitudes 
with the usual scalings dt ~ 0(A 2 ) and V ~ O(Ai), pro- 
vided that the coefficient of the quadratic resonant term 
is small. After rescaling they can be written as: 
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The equations for the other two amplitudes are obtained 
by cyclic permutation of the indices. The broken chiral 
symmetry manifests itself by the cross-coupling coeffi- 
cients not being equal. The coefficient 7 is, therefore, a 
measure of the rotation. 

A study of the steady hexagons with a wavenumber 
k slightly different from the critical wavenumber (k„ = 
k-n + q n , |<Jn| = l) an d their sideband instabilities has 
been undertaken in jl6| . We focus here on the secondary 
Hopf bifurcation that appears for larger values of the con- 
trol parameter fj, = (i c = (2 + v)/{v - l) 2 + q 2 RQ. It 
gives rise to what we call oscillating hexagons, in which 
the three amplitudes of the hexagonal pattern oscillate 
with frequency to — 2^/2>^j/{v — l) 2 and a phase shift of 
27r/3 among them. As \x is increased further, eventually 
a point // = fi^et is reached at which the branch of os- 
cillating hexagons ends on the branch corresponding to 
an unstable mixed-mode solution in a g lobal bifurcation 
involving a heteroclinic connection |l7]-[l^]. Above this 
point the only stable solution is the roll solution. When 
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I7I > v — 1 the rolls are never stable and the limit cycle 
persists for arbitrarily large values of /i. In the absence 
of the quadratic term in Eq. (|l|) this condition corre- 
sponds to the Kiippers-Lortz instability of rolls [^0| . Far 
above the Hopf bifurcation the periodic orbit is expected 
to become anharmonic, and may somewhat resamble the 
state encountered in the Kiippers-Lortz regime of rotat- 
ing Boussinesq Rayleigh-Benard convection |^l],^2| . 

To study the stability of the oscillating hexagons close 
to the Hopf bifurcation point, the amplitudes A n are ex- 
panded as: 
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The amplitude of the steady hexagons at the bifurcation 
point, R — — 1), is independent of p and q. Since we 
are considering a secondary bifurcation the amplitude Ti 
will couple to the phases of the hexagonal pattern [p3j . 
From the phases 4> n of each of the modes it is possible to 
construct a phase vector tf> = (<p x , <f> y ) = (—(f>2 — 4>3, {<p2 — 
(f>3)/V3), whose components are related to translations 
in the x- and y-directions |p4]-p6||. Eliminating the fast 
variables, we arrive at order at an equation for the 
amplitude of the oscillation H, coupled to the phases of 
the underlying hexagonal pattern, 
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where e z is the unit vector in the direction perpendicular 
to the plane, e = fj, — /j c , u) = 2\/3"/R 2 and 

v = 3R(1 + 2R), 
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It is worth pointing out that the phase-amplitude equa- 
tions can be deduced by means of symmetry argu- 
ments alone and are, therefore, generic to this order in 
e. In fact, they could be derived directly from the fluid 
equations without the use of the Ginzburg-Landau equa- 
tions (|l|). Thus, keeping higher order terms in (|l]) would 



change the values of the coefficients in , but not their 
form. 

Central to the results presented in this Letter is the 
observation that for hexagons with rotation as described 
by ([l]) the phase-amplitude equations (0,0) decouple at 
the bandcenter, i.e. for q = 0. In this case they reduce 
to the usual CGLE for the amplitude of the oscillation. 
Scaling the time, space and amplitude, Eq. (|J) can be 
written in the more usual form |q| : 

d t H = H + (1 + ih)V 2 H -{b 3 ~i sign(uj))H\H\ 2 (5) 

where b\ = £i/£ r = and 63 is given by: 
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As uj is increased by increasing the rotation rate 7, 
63 reaches a maximum 6g nax = 6|j nax (jy) at ^ max = 



y/2(v — l) 2 (^ + l)/3/(^ + 3). For any value of v the co- 
efficient 63 will move across the range < 63 < bf & * 
as the rotation rate is varied. Note that the depen- 
dence of 6™ ax on v is only weak and is limited to the 
range 0.354 < bf™ < 0.375, as indicated by the dotted 
lines in Fig. |l|. Thus, independent of v 1 which depends 
on the physical properties of the system, the oscillating 
hexagons are always in an intermittent regime in which 
stable plane waves and defect chaos coexist (see Fig. [I]). 
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FIG. 1. Phase diagram for the CGLE (after [6]). Defect 
chaos is found to the left of line T. Above the Benjamin-Fair 
(BF) line, all plane waves are longwave unstable, resulting in 
a state of phase chaos, which is bistable with defect chaos or 
a frozen vortex state up to the line L. The lines Si and 5*2 
represent the convective and absolute stability limits of plane 
waves emitted by spirals. The vertical dotted lines represent 
the limits for the range of &3 lax (cf. (^)). As q is increased 61 
decreases. 

It should be noted that, although we are considering 
a system with broken chiral symmetry, to this order in 
the perturbation expansion we obtain the usual CGLE, 



2 



which is chirally symmetric. Going to higher order we 
would expect terms like e z ■ (VH x VH*)H to appear 
in Eq. (^), breaking the symmetry between spirals with 
opposite topological charge fl27[. 



the influence of this coupling in the defect chaotic regime 
we measure the density of defects for Eqs. ([!]) and Eqs. 
([|[|) as a function of q (accordingly bi) and compare it 
with the results for the CGLE (§). 




Recon- 



with 



FIG. 2. 

stmction of the hexagon pattern ip — y^ 3 _ ] Ajt 
k c — 30ir/L x , obtained by simulating numerically Eqs. (|l|) in 
a box of length L x = 50, L y = 100 /VS with 128 x 128 modes, 
for fj, — 4.6, v = 2, 7 = 0.5 and q = 0. The contour lines 
are taken at ip = —0.7. The time difference between the two 
snapshots is half a period of oscillation. 

We investigate the defect chaos regime by numerical 
simulations of Eq. ([j]) using a pseudospectral method 
with a 4th-order Runge-Kutta/integrating factor time- 
stepping scheme and periodic boundary conditions. To 
allow for regular hexagonal patterns we take a rectan- 
gular box of aspect ratio L x /L y = \[2>j1. Fig. || shows 
a picture of the hexagonal pattern in this regime. It 
consists of patches of slightly roll-like hexagons, whose 
preferred direction oscillates on a fast time scale. The 
patches change shape and size on a slow time scale. The 
regions where almost perfect hexagons can be observed 
(e.g. on the bottom left part of the figure) correspond 
to the zeros (defects) of the oscillation amplitude. In 
order to extract this complex amplitude Ti. from the am- 
plitudes A„ we use that close to onset He iut + H* e~ iut ~ 
\A n \ — (Y^j—i \Aj\)/3. The amplitude TL is obtained mul- 



tiplying the former expression by 



and taking the 



average over each period. A snapshot of the magnitude 
\TL\ as obtained from Fig. || is given in Fig. ||a, while in 
Fig. [|d we show the corresponding lines Ke(TL) = and 
Im(7Y) = 0, whose intersection points correspond to the 
defects of TL. 

Away from the bandcenter, q ^ 0, b± becomes non-zero 
and is given by: 
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The coefficient 63 remains unchanged. Depending on the 
value of q the system can therefore cross the line T (see 
Fig. [l]) , beyond which defect chaos is no longer an attrac- 
tor within the CGLE (||). But, for q 7^ 0, the oscillation 
amplitude TL is coupled to the phase (p. In order to study 





FIG. 3. Snapshots corresponding to Fig. |^ of the a) modu- 
lus of TL, b) lines Re(H) = (solid) and lm(H) = (dashed). 



We begin with a perfect steady hexagonal pattern with 
wavenumber q (A n — i?e lq "' x , |q„| = q, in Eq. (Q), 
TL = 0, 4> = in Eqs. (|,|), H = in Eq. (|)) and 
add noise of zero mean. If the system is large enough 
the resulting oscillating state is not homogeneous but 
ends up in a persistent chaotic state. We also checked 
numerically that a perfect oscillating state in Eq. (Q), 
A n = (i?+We iw * + c.c.)e iqn ' x , is linearly stable, although 
sufficiently strong perturbations can destabilize it, giving 
rise to defect chaos. In the simulations we use 64 x 64 
modes, in a box of length L x = 200 for Eqs. (@) and (|J|), 
and L = 50 for the CGLE, which roughly corresponds to 
the former after rescaling. As the scaling depends on q 
we take the same box size L Xl L y in all the simulations 
of Eqs. (Q) and (||,|]) and then scale the density of de- 
fects appropriately to compare with the results from the 
CGLE. After a transient time, the number of defects is 
measured for a sufficiently long time so the system is 
statistically stationary. The results corresponding to the 
CGLE are the average of three independent runs. 

As can be seen in Fig. [| the simulations of (|l|) , (^Ji)) 
and (||) agree very well up to the line S2 , where the spiral 
defects go from being absolutely to convectively unstable. 
Since for q ^ the amplitude Ji is no longer decoupled 
from the phases, the agreement indicates that the phase 
does not have a strong influence on the dynamics in this 
regime. Beyond line S2, after a chaotic transient, the 
system settles down in a frozen vortex state in which the 
density of defects remains unchanged for long periods of 
time. Since this final state is history dependent a quan- 
titative agreement is harder to achieve. In particular it 
would require a more extensive sampling of initial con- 
ditions. However, the agreement remains good. Beyond 
line T still a small number of defects remain, probably 
because our simulations are not long enough given the 
weak interaction of defects. If and to what extent the 
coupling with the phase affects the lines S2 and T is not 
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clear from these simulations, but the qualitative picture 
remains the same as in the CGLE. Note also that, for 
these parameters, additional long-wave instabilities oc- 
cur beyond line T |28| . 
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FIG. 4. Density of defects n as a function of the wavenum- 
ber of the underlying pattern q for 63 = 0.355 (y = 2, 7 = 0.5, 
H = Uc + 0.1). The circles correspond to the simulations of 
Eq. (|5j), while the squares have been obtained simulating Eqs. 
(^,^). The stars are the results obtained with the amplitude 
equations (jjj) . Between lines 5*2 and T the system ends up in 
a frozen state. 

In conclusion, we have shown that within weakly non- 
linear theory hexagons arising in rotating convection 
are a good candidate to investigate two-dimensional de- 
fect chaos. Due to the rotation the hexagons undergo 
a Hopf bifurcation to oscillating hexagons. Although 
these are usually linearly stable at the band center, fi- 
nite perturbations that excite defects lead to persistent 
spatio-temporal chaos. In this regime the single complex 
Ginzburg-Landau equation describes the dynamics quan- 
titatively. Further away from the band center we find a 
transition to a frozen vortex state, which is also in agree- 
ment with the CGLE. The decoupling of the phases and 
the amplitude that occurs at the band center will be mod- 
ified if higher-order corrections are taken into account in 
the amplitude equations (1). Simulations in which non- 
linear gradient terms in (1) are included suggest, how- 
ever, that the qualitative picture remains the same [p8[ . 
It is worth noting that the dynamics discussed in this 
paper (in particular, the spatio-temporally chaotic state) 
have been obtained for values of the rotation rate be- 
low the Kiippers-Lortz instability. Therefore we expect 
that this genuinely new regime of rotating convection is 
accessible in currently available experiments. We hope 
that these results will trigger new experiments, which 
will contribute to a better understanding of the mecha- 
nisms of transition to spatio-temporal chaos. 
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